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This is the second of two companion papers on computing the self-force in a radiation gauge; 
more precisely, the method uses a radiation gauge for the radiative part of the metric perturbation, 
together with an arbitrarily chosen gauge for the parts of the perturbation associated with changes 
in black-hole mass and spin and with a shift in the center of mass. In a test of the method delineated 
' in the first paper, we compute the conservative part of the self-force for a particle in circular orbit 

around a Schwarzschild black hole. The gauge vector relating our radiation gauge to a Lorenz 
gauge is helically symmetric, implying that the quantity h a pu a u^ must have the same value for 
our radiation gauge as for a Lorenz gauge; and we confirm this numerically to one part in 10 14 . 
C ' As outlined in the first paper, the perturbed metric is constructed from a Hertz potential that 

is in term obtained algebraically from the the retarded perturbed spin-2 Weyl scalar, ip r $ . We 
use a mode-sum renormalization and find the renormalization coefficients by matching a series in 
^-l L = I + 1/2 to the large-L behavior of the expression for the self- force in terms of the retarded field 

Kxh we similarly find the leading renormalization coefficients of hapu 01 ^ and the related change in 



5-H 



the angular velocity of the particle due to its self-force. We show numerically that the singular part 
of the self- force has the form = (V a p _1 ), the part of V Q p _1 that is axisymmetric about a radial 
line through the particle. This differs only by a constant from its form for a Lorenz gauge. It is 
because we do not use a radiation gauge to describe the change in black-hole mass that the singular 
' part of the self-force has no singularity along a radial line through the particle and, at least in this 

\& ' example, is spherically symmetric to subleading order in p. 

' 

00 ' PACS numbers: 04.30.Db, 04.25.Nx, 04.70.Bw 

o\ 
o 

5 ' I- INTRODUCTION 



We present here a first self-force computation in a radiation gauge, following the method outlined in a companion 
paper |l[ (henceforth Paper I) . The computation has been done previously by Barack and Sago Q and by Detweiler 
[3|, and gauge invariant quantities associated with the conservative part of the self-force are compared in their joint 
paper A radiation-gauge approach has the advantage that one can use the Teukolsky equation to compute the 
perturbed metric and the self-force for orbits in a Kerr background, and the present paper serves as a test of methods 
described in Paper I for that problem. 

The use of a radiationgauge for the self-force problem has been delayed in part because the MiSaTaQuWa renormal- 
ization prescription 0, [6|] was developed for a Lorenz gauge and in part because, in a radiation gauge, the linearized 
metric of a point-particle is singular along a ray through the particle. One can avoid a singularity of this kind by 
restricting the use of a radiation gauge to the part of the perturbation determined by the gauge-invariant Weyl scalar 
ipo (or ipi). The part of the metric perturbation that describes the change in mass and angular momentum of the 
spacetime can then be added in any convenient gauge. The perturbed metric obtained from from ip is constructed as 
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a sum of angular and time harmonics, defined for r > tq and for r < ro, with ro the Schwarzschild radial coordinate 
of the particle. 

Although the £ > 2 part of the metric perturbation can be computed more simply in a radiation gauge, an analytic 
computation of the singular field that is to be subtracted is significantly more difficult. We avoid the difficulty by 
replacing the analytic computation by a numerical determination of the renormalization coefficients that are subtracted 
in a mode-sum renormalization of the self-force. As a result, the efficacy of the method depends on the numerical 
accuracy with which these coefficients can be computed. We describe the numerical methods used and report tests of 
their accuracy. 

As noted in Paper I, recent work by Gralla 0, following an earlier derivation of the self-force equations by Gralla 
and Wald 0, shows that the first order correction to the geodesic equation, 

u p V p u a = a rcnQ , (f) 

is obtained from 

« retQ = -(9 aS - u a u s ) (vpK$ ~ ^5hfA u^u\ (2) 

by an angle average in locally inertial coordinates, over a sphere of geodesic radius p about the particle: 

a"" " = lim f a rct "dn. (3) 
Js p 

The equation holds in any gauge for which the leading part of h T £g is 0{p~ l ) and has even parity. Paper I showed that 
the even-parity condition was satisfied in a radiation gauge. Taking the angle average is equivalent to subtracting a 
field a s (the singular part of the acceleration) 

„rena ret a sa / a \ 

a = a — a , (4) 

if a sa satisfies the conditions (i) The limiting angle average of a s , defined as in Eq. ([3]), vanishes; and (ii) a TCta — a sa 
is continuous at the particle. 

The numerical determination of a sa from a rctQ shows, with an accuracy close to machine precision, that a sa is 
proportional to (V a p _1 ) and hence is spherically symmetric. The result also implies that the singular field can be 
identified with its leading and subleading terms in its mode-sum expression as a power series in L = £+ 1 /2. It remains 
an open question whether this unexpectedly simple behavior of the singular part of a a holds for generic orbits or for 
a Kerr background. 

The plan of the paper is as follows. The perturbed metric is constructed from the Weyl scalar ipQ Ct , computed as a 
sum over spin- weighted spherical harmonics, and Sec. HI] details the numerical method used to compute the retarded 
radial and angular functions that comprise the sum. In particular, radial integrations using the Teukolsky radial 
equation and the Sasaki-Nakamura form are used and compared. Sec. IIIII describes the computation of the retarded 
metric and the retarded expression for the self force from the values of ipo for each harmonic. The conservative part of 
the self force has only a radial component, f r , and it is renormalized by matching a power series in L to the sequence 
of contributions f r [£] from successive angular harmonics. We find that the singular field obtained in this way is the 
angular decomposition of (V"/)" 1 ), with p the geodesic distance orthogonal to the particle's trajectory. 

Associated with the perturbed metric of a particle in circular orbit is a set of related quantities that are invariant 
under helically symmetric gauge transformations. The Sago et al. comparison paper Q tabulated values of one of 
these quantities, and we add the values from our radiation-gauge computation to their comparison table. 



II. COMPUTATION OF ip T et 

Formalism 

We begin with a brief review of the formalism used in Paper I. We consider a particle of mass m in circular orbit 
about a Schwarzschild black hole of mass M. We use Schwarzschild coordinates, with the particle at radial coordinate 
ro in the 9 = 7r/2-plane. The particle's four- velocity is 



u a = u t (t a + ncj) a ), 



(5) 
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with t a and <p a timelike and rotational Killing vectors, il = y/M/r$, and it* = 1/yT — ZM/tq, With a <5-function is 
normalized by / S A (x, z)^ [ /\g\d 4 x = 1, its stress-energy tensor is given by 



T 



a/3 



mu a U P / 5 A {x a -z a { T ))dT 



u a u p 



m- 



u l (-g) 



m 5{r-r Q )5{cos6)8{4>-Slt) 



(6) 



where a change of coordinates from r to t in the integral for stress-energy tensor is used to obtain the second equality. 
From the addition theorem for spin-weighted spherical harmonics, we have 



T afi = ^2 



Lm 



mu a u^ 
r 2 u l 



5{r - t q ) s Y em (0, </>) a Y tm (-K/2, Sit). 



The Kinnerslcy tetrad vectors have components 

(I") = (l//(r), 1,0,0), K) = i(l,-/(r),0,0), (m") = -±= 

£ v 2r 



(O,O,l,i/sin0), 



(7) 



(8) 



where f(r) := 1 — 2M/r; and we denote by D, A, and 6 the derivative operators along the tetrad vectors I, n and m, 
respectively. The nonvanishing spin-coefficients associated with this tetrad arc 



q = — 1/r, P = —a 



cot 6< 
2V2r 



A 

2^' 



and 7 



M 
2^' 



(9) 



where A = r 2 — 2M r. 

The pertubed Weyl scalar tpo = —C a p~ l sl a mPl 1 m s satisfies the Bardeen-Press equation (the Teukolsky equation for 
o = 0), 



A * 



f„2 



Mr 



r)d,-A- 



8 







dr V dr 



i'o = 4irr 2 T, 



where 



T := -2(8 - 2P)8Tu + 4(D - 4g)(« - 2/3)T 13 - 2(D - 5g)(D - g)T 33 
= t ( °) +T {1) + T (2) . 

The superscripts indicate the maximum number of radial derivatives in each term. Explicitly, 



T(0) = _y-m^ J(r _ ro)[ ^_ 1 )t(i + 1 )(i + 2 )]V\Y tfn (6,<f>)Y tm (?-, 



nt 



r (D = ^ 



2z<5'(r-r ) + — — <5(r-r ) 
Jo 



[(* - 1)(£ + 2)] 1 ' 2 2 Y hn {0, faYtn (|, fit 



(10) 

(11) 

(12a) 
(12b) 



T< 2 ) = ^mfl 



V 



S"(r - r ) + - - 



6 2i?nO\ , 



/o ; 



*'(r - r ) 



i 2 2 6imft 2imMfi\ 



/o n)/o 



<5(r - ro) 



■ 2 Yi m (9, cp) 2 



Yt, 



-,nt 



The Weyl scalar has the harmonic decomposition 

•0o = e ~ iUmtR o^ m ( r ) 2Y em (e, 0), cj m := mSl. 



(12c) 



(13) 



(in 



Its radial functions Roeu m and the corresponding radial functions Ruu m of p 4- 04 satisfy (^ and uj m subscripts 
suppressed) 



Ai? ' + 6(r- A/)i? + 
A< - 2(r - M)i£ + 



w 2 r 4 4iwr 2 (r - 3M) 
~A~ + A 
wV 4iwr 2 (r - 3M) 



- (*-2)(* + 3) 
-(*-!)(* + 2) 



i? = 0, 
R 4 = 0, 



(14) 
(15) 
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where the prime denotes a derivative with respect to the radial coordinate r. Solutions to the above equations are 
related by 



in this relation, however, Rq and i?4 differ by a relative normalization from the radial functions of components ipo 
and ip4 °f the same vacuum Weyl tensor. 

We denote by and R^ solutions to Eq. (|14[) that are regular at the horizon and at infinity, respectively, and 

will write (') := — . As shown in Sect. V of paper I, the solution to Eq. (flUl) is given by 
dr 

= 4 0) +4 1) +4 2 \ (17) 

where the three terms, corresponding to the three source terms of Eq. ()12l) . have the form 

4 0) = 4W^^A fm [(€-lK(£+l)(^ + 2)] 1 / 2 i? H (r < )i? oo (r > ) 2 r, m (0,0)r, m (|,m), (18a) 

Im 

£m 

{[imfirg + 2r ]R li (r < )R oo (r > ) + A [ J Ri i (r )i? oo (r)% - r ) +i? H M#>o)%o -r)]}, (18b) 

4 2) = -47rmnV^^ m2 y /m (e,^) 2 yi m (|,m) x 

[30r^ - 80A/rjj + 48A/ 2 r 2 - m 2 r2 2 r£ - 2A 2 , - 24A r (r - M) + 6im^(r - A/)]i? H (r < )i? 00 (r>) 
+2(6rg - 2QM4 + l6M 2 r 3 - 3r A 2 + irr&Aort)[Ik(ro)Roo{r)0(r - r ) + iC(r o )Jfe(r)0(ro - r)] 
+r 2 A 2 [<(r o )i? oo (r)0(r - r ) + iC(r o )J2H(r)0(r o - r) + W[J2 H (r), i?oo(r)]<5(r - r )]|. (18c) 

Here the Wronskian-related quantity, 

Alm := ^{RkR'^-RooR'kY (19) 

is a constant, independent of r. 
Numerical method 

To compute tpQ Bt , we use a 7th order Rungc-Kutta routine to integrate the radial Teukolsky equation, matching the 
radial function to a power series expansion near the horizon and infinity. Because our renormalization method relies 
on numerical extraction of the renormalization coefficients, we work to high numerical precision and check the results 
by comparing independent codes based on the Teukolsky form of the radial equation (fl"4"|) and on the Sasaki-Nakamura 
form. We find the spin-weighted spherical harmonics s Yg m (^ ,0) to similarly high precision. 

Because one of solutions for Rq diverges at the horizon, we integrate R4 from the horizon and Rq from infinity. We 
match near the horizon and at large r, respectively, to the series expansions 

i? 4 H = e^J2*n(jj-2)\ (20) 

n 

Rooo = e lUJr ( r /M)n+5 ' ( 21 ) 

where r* = r + 2M log(r/2A/ — 1). The expansion coefficients are found from a 3-term recurrence relation, 

-2iwM {n - 5)a„_ 2 + (I 2 + I - 6 + 5n - n 2 - 8incoM + 20iwA'/)a„_i 
fl " = 2n{n - 2 + AiujM) ' ( } 

((6 + 8n + 2n 2 )6 n _ 2 + (f + £ - 2 - 3n - n 2 )&„_i) , (23) 



2nuiM 
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obtained by substitution in Eqs. (|T4|) and (|T5l) . The integrations from the horizon and from infinity yield the two 
independent solutions to the homogenous Teukolsky equation that we have labeled i?n and Roq. By using the outgoing 
wave solution at the event horizon as an initial condition for , we get the correct boundary condition (ingoing at 
horizon) for s = +2 radial solution on using Eq. (fit)]) i.e. 



R "-^JT- ( 24 ) 



The numeric results are compared with the numerical solution obtained from Sasaki-Nakamura equation for con- 
sistency. In this case, we integrate 



d*X t , 



dr* 

where 



= U(r)Xi >mfi> , (25) 



TT , , 12M 2 -2Mr(£ 2 + e + 3)+r 2 (£ 2 + £-r 2 uj 2 ) 

U(r) = 3 . 26) 

The function Xe tm ^ is related to the homogenous solution of the s = —2 radial, Teukolsky equation by 

2A(r-3M + w-V) , /Z(Z + 1)A 6A/A 2ru{3iM - ir + r 2 u)\ 

4,Re,m,u — X t + - Xi mu! , (27) 

rjr ' ' \ rjr r/r z r\ J 

where r) = (I — 1)1(1 + + 2) — 12iMoj. We use a 7th order Runge-Kutta routine to integrate Eq. ([25|) with the 
boundary conditions 

X£ m>u = e^E^G^" 2 )"' (28) 

n 

= ^*E(^w (29) 

where c„ = d n = for n < 0. The values of c„ and d n are calculated from the following recurrence relations: 

i(n-3)Mw e 2 +£- (n-2)(n-3 + 12iMw) . , 

C "-3 H T~? , A ., T \ c n-2, (30) 



2n(n + 4iMuj) An{n + AiMu 

£ 2 + i - 2n 2 + bn - 6 - 12i(n - l)Mu 



2n(n + 4iMw) 



" Cn— 1 j 



. (n-3)(n+l) J . (l + n )(l-n + l) , 
d " = * 2»M, d - 2 - * 2^ ^ (31) 

The outgoing (ingoing) boundary condition at infinity (event horizon) is used in combination with Eq. (|16p which 
relates the s = 2 and s = — 2 radial solutions. 

The accuracy of the homogenous solutions is monitored by the constancy of Aj = A 3 W(Rn, i?oo)- Table U shows 
the fractional standard deviation (f.s.d.), or the departure of A 3 !^ from its average value. A comparison of the 
s = ±2 radial functions given by the Teukolsky equation and the Sasaki-Nakamura equation (Eqs. (fl4| and (fTS"])) arc 
included. Measured in this way, direct integration of the Teukolsky form yields 13-digit accuracy. Table II reports a 
second check of numerical accuracy by directly comparing the radial functions computed from the two forms of the 
radial equation. 
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I 


Mlo 


f.s.d.(T-eqn) 


f.s.d.(SN-eqn) 


3 


0.001 


4.2 


X 


1Q -14 


2.8 


X 


1Q -14 


2 


0.136083 


2.8 


X 


1Q -14 


9.7 


X 


IO" 13 


5 


0.000353553 


4.4 


X 


io- 14 


2.7 


X 


10 -14 


5 


0.340207 


3.5 


X 


io- 14 


1.3 


X 


io- 12 


10 


0.00108866 


5.6 


X 


1Q -14 


4.3 


X 


10 -14 


9 


0.612372 


4.3 


X 


io- 14 


1.0 


X 


io- 12 


14 


0.000715542 


6.2 


X 


io- 14 


5.5 


X 


io- 14 


15 


1.02062 


5.1 


X 


10 -14 


1.1 


X 


io- 12 


20 


0.000929429 


7.8 


X 


1Q -14 


7.4 


X 


10 -14 


21 


1.42887 


5.9 


X 


1Q -14 


1.4 


X 


IO" 12 


25 


0.000707107 


9.2 


X 


1Q -14 


9.7 


X 


10 -14 


25 


1.70103 


7.2 


X 


io- 14 


1.1 


X 


io- 12 


30 


0.002 


1.0 


X 


IO" 13 


1.2 


X 


IO" 13 


31 


2.10928 


8.0 


X 


IO" 14 


1.6 


X 


10 -12 


35 


0.000544331 


1.2 


X 


IO" 13 


1.2 


X 


IO" 13 


34 


2.31341 


8.9 


X 


1Q -14 


1.1 


X 


IO" 12 


40 


0.000988212 


1.3 


X 


IO" 13 


1.5 


X 


IO" 13 


40 


2.72166 


9.6 


X 


io- 14 


1.6 


X 


IO" 12 


45 


0.000603682 


1.5 


X 


IO" 13 


1.7 


X 


IO" 13 


46 


3.1299 


1.1 


X 


io- 13 


1.3 


X 


io- 12 


50 


0.000353553 


1.6 


X 


IO" 13 


1.8 


X 


IO" 13 


50 


3.40207 


1.2 


X 


IO" 13 


1.5 


X 


IO" 12 


55 


0.00067466 


1.8 


X 


io- 13 


2.0 


X 


io- 13 


55 


3.74228 


1.3 


X 


io- 13 


1.1 


X 


IO" 12 


60 


0.001 


2.0 


X 


IO" 13 


2.3 


X 


IO" 13 


61 


4.15052 


1.3 


X 


io- 13 


1.5 


X 


nr 12 


65 


0.000544331 


2.2 


X 


i<r ia 


2.5 


X 


io- 13 


64 


4.35465 


1.4 


X 


IO" 13 


1.5 


X 


IO" 12 


71 


0.000494106 


2.3 


X 


io- 13 


2.8 


X 


io- 13 


70 


4.7629 


1.7 


X 


IO" 13 


1.4 


X 


IO" 12 


75 


0.000637528 


2.5 


X 


IO" 13 


3.1 


X 


IO" 13 


74 


5.03506 


1.8 


X 


io- 13 


1.5 


X 


io- 12 


80 


0.000471818 


2.8 


X 


IO" 13 


3.5 


X 


IO" 13 


80 


5.44331 


1.9 


X 


10" 13 


1.5 


X 


IO" 12 


85 


0.000707107 


3.0 


X 


nr 13 


3.9 


X 


io- 13 


85 


5.78352 


2.1 


X 


nr 13 


1.0 


X 


IO" 12 



TABLE I: Accuracy of radial integration measured by fractional standard deviation (f.s.d) of A 3 !^ from its average value, for 
specified I and u). The third and fourth columns list the f.s.d. obtained by integrating the Teukolsky and Sasaki-Nakamura 
forms of the radial equation, respectively. The frequencies are chosen to lie between the maximum and minimum of mil, 
between the value Mu> = £/6 s/2 at the ISCO and Mu ~ 1/150 3/2 , at r ~ 150M, m = 1. The f.s.d. in the region r = 6M to 
20M is usually twice the average shown above. 



To calculate s Yf m to high precision, we used the following analytical forms of spin- weighted harmonics at 6 = tt/2. 
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r /M 


£ 


m 


fractional difference 


100 


2 


1 


5.9xl0~ 13 


10 


2 


2 


4.6xl0~ 15 


6 


2 


1 


1.9xl0~ 13 


80 


75 


75 


1.7xl0~ 15 


80 


20 


1 


3.2xl0~ 15 


13 


5 


4 


1.6xl0~ 14 


10 


3 


1 


2.3xl0~ 14 


8 


15 


14 


5.0xl0~ 15 


8 


12 


6 


1.6xl0~ 14 


70 


3 


2 


1.8xl0 -15 


6 


4 


1 


4.1xl0~ 13 


50 


4 


1 


l.OxlO" 13 


6 


10 


9 


4.7xl0" 13 


70 


9 


8 


1.7xl0~ 15 


6 


25 


1 


4.0xl0" 16 


7 


25 


25 


2.4xl0~ 15 


75 


19 


18 


5.0xl0~ 15 


6 


20 


20 


2.0xl0~ 14 


72 


85 


78 


1.5xl0~ 15 


6 


85 


85 


2.4xl0~ 14 


7 


85 


5 


4.2xl0~ 14 


60 


75 


1 


4.0 x 10" 16 


10 


75 


71 


2.2 x 10~ 15 


15 


50 


38 


1.4 x 10" 15 


6 


40 


25 


7.2 x 10" 15 


6 


25 


15 


7.1 x 10~ 15 


30 


25 


1 


6.0 x 10~ 16 


20 


25 


25 


7.8 x 10" 15 


10 


2 


1 


3.9 x 10~ 14 



TABLE II: For each listed value of £, m and ro, we give the fractional difference between the radial functions at ro, obtained 

i/2 3/2 

by integrating the Teukolsky and the Sasaki-Nakamura forms of the radial equation, with frequency lu = mQ, = mM 1 r . 



Introducing the symbol e^ m := 



1 £ + m even 
0, I + m odd , 



we can write 



V 2' ' y ' V^(£-m)U(£ + m)\\ 



(£+m)/2 (2£+l)(£-m)l(£ + m)\ 



met m 



2^ m (-,0) = (- 



(£+m)/2 (2£+l)(£-m)\{£ + m)\ 
V Att(£ -!)£{£ +!){£ + 2) 



{£ - m)\\(i + to)!! (£-m- + to - 1)!! 
[2m 2 - £(£ + 1)] e e , m 2ime (tm+1 



(£ - m)\\{£ + to)!! (I - to - !)!!(£ + m - 1)!! 



(32) 

(33) 
.(34) 



In this way, we obtain values of ?pQ et with accuracy of 1 part in 10 



III. METRIC PERTUBATION AND SELF-FORCE 



We use the outgoing radiation gauge (ORG) satisfying the conditions 

h af3 n a = 0, h = 0. 



(35) 



8 



Wc find the retarded perturbed metric from a Hertz potential <J ,ret satisfying 

8$j et = g 4 * rct + 12Md t $ ie \ 
whose algebraic solution for each angular harmonic is given by 

c (-l) m (i + 2){l + !)£{£ - l)$ t ,-m + UimMfl^m 



[(£ + 2){£ + 1)£(£ - l)] 2 + lUm 2 M 2 n 2 

±\ JD —imVlt 



where = J2 Lm * im (r) 2 Y irn (6, 0) e -"» ot and Vo = Et, m lM»0 2 Y t . 
The ORG form of the metric pertubation h a p is 



h a t 



r 4 { n a n p {8 + 2/3) (5 + 4/3) + m a mp(A + 5/i - 2 7 )(A + fi- 4 7 ) 



-ri( a mp) 



[{8 + 4/3)(A + ^ - 4 7 ) + (A + 4 M - 4 7 )(5 + 4/3)] }tf + c.c. 



(36) 



(37) 



(38) 



where the spin coefficients are given by Eq. ([9]). After writing d in terms of S, the tetrad components take the form 



/in 
h 3 3 
his 



d 2 - 2/9 t d r + / 2 d 2 3(r-M)_ , /(3r-2M) 



2r 2 



2r 2 



2M 2 



2\/2 

Finally, the self-acceleration, 



fit - fd r - - 



-(9 



a. 5 



u A ) V h 



rcn 

-yS 



(39) 
(40) 
(41) 

(42) 



written in terms of these tetrad components, is given by 



1 M 



— f D + -A + -fi(S-3) -- — 



16"" ' 8 ' 8""'" 2r 2 

+ / -— £>--_ A + -— )h 33 

[\8r 4r / 2 r 2 J 



1 M 
4%/2 r 2 lC 



2\/2 



^ ^13 + C.C. 



(43) 



The harmonic decomposition, VP = $ imir) 2^ 



of the Hertz potential gives a corresponding decom- 



position of a r . The contribution to a r from the harmonic ^£ m , however, is not a single angular harmonic because 
the particle's velocity u a involves the Killing vector <fi a , an axial £ = 1 vector field on the 2-sphere; and u a occurs 
quadratically in Eq. (|42|) for a r . Instead, the terms in Eq. (|43f include terms from hu and 8/113 proportional to 
Yi m ; terms from Qhu and /113 proportional to iY^ m ; and terms from 8/113 and /133 proportional to 2^?™. A virtue of 
our numerical rcnormalization procedure is that we need not rewrite these latter terms as sums of spin-weight zero 
harmonics. 

We organize the computation by writing the right side of Eq. (|43[) as a sum of six terms for which the magnitude 
I s I of the spin- weight increases from one red square bracket to the next: 



1 , „ 3 . 1 M 
—f Q D + -A-- — 

16 8 2 r„ 



hu 



1 



4\/2 



dhi 



1 , 



^(«*) a n/o(roA-i)(fc ls -ft 14 ) 



.4^2 

6 



( U *) 2 ^(^13 + 8/U4) + (U*) 2 ^ 



8 4/o 



1 

2^ 



(/l33 + ft.44) 



(44) 
(45) 



1 = 1 



9 



The subscript i (i = 1, ... ,6) in the symbol a\ refers to the location of the term on the right side of Eq (|44[) . We 
denote by a r [£] the contribution to a r from the restriction of VP to the £-subspace, writing 



(46) 



The harmonics a[[£] have finite limits as r 



We evaluate a r at a particle position P with coordinates t = 0, 



6 = it/2, and 
forms 



0, where the value of Mm is real. The harmonics then have at P (for either choice of limit) the 



al\£] 
a r S) 



■[{£-!)£{£ + !)(£ + 2)] 1 '*(u t ffffi 



2d t + f [d r + - + 4 



M 



- 1K(£ + 1)(£ + 2)} 1 / 2 (u t ) 2 Mr f £ 

m 

—£(£ + !)[(£- !){£ + 2)} 1 / 2 (u t ) 2 fy Q nY,^im)(iY em + -Mm), 



Rc(^ £m ) 



Ytm, 



[{£ - !){£ + 2)] 1 / V) 2 / r 4 ft ( % ~ 2f d t d r + f 2 d 2 



— (1 + /o)5 t H 



3-2- 



M 

l 



9 r 



f£m, 



+ -|(l + 2/ )J J2 Im ^™ft 

^ / TO 

1)(* + 2)( !t i ) 2 / Mr E( mfflm (*^) ~ /od,.Rc(*£ m ) - 2/rRc(tf fm )) 2 Yfa 



dfd r - 2f d t 8 2 + f 2 df, + -8 2 - 2(9 - 13M/r )d t d r + 12^(1 - AI/r )d 2 

ro r 



-|(5 - 4M/r )d t + 4(17 - 32M/r + 8A/ 2 /r 2 )c> r + if (1 - M 2 /r 2 ) 
rn rn Tn 



(47a) 
(47b) 
(47c) 
(47d) 

(47e) 



E M^ira) Mtm- (47f) 



The mode-sum renormalization of the self-acceleration is given by 



n ™ E( c 



(48) 



f=0 



with a s r [^] the singular part of a rct r [£] . As described in paper I, one can identify a sa [£] with the leading and subleading 
terms in L in the large- L expansion of a retQ . This expansion has the (direction-dependent) form 



oo 



fc=l 



(49) 



where cn falls off at large £ faster than any power of £, and the superscript ± refers to the limit r 7q . Because 
we find that the singular field can be identified with its leading and subleading terms, the remaining part of the 
power-series expansion in L must sum to zero, and we reorder the higher-order terms in the power series, replacing L k 
by a sequence of polynomials Pzk {£) that individually sum to zero and whose leading term is L k (see Eq. ([HI), below). 
We find the singular field numerically by matching a power-series of this form, truncated at a finite value fc ma x of k, 
to the sequence a ret [£], computed using Eqs. (f4"6")l and (f4"7| . Since the £-mode expansion of the singular field agrees 
with the largc-£ expansion of the retarded field, one can extract the regularization parameters by the above method. 
Though we find the singular field ipf, to order 0(e _1 )in the previous paper, one can see that it takes a heavy amount 
of analytic work to calculate each of the af r [£] which requires a careful analysis of more than 200 terms! Hence, we 
employ this numerical matching method to extract the renormalization coefficients in the radiation gauge. 

A striking feature of the numerically determined singular field is that the first three terms in the L expansion of the 
conservative part f r of the self- force coincide with the expression for a B r obtained from the power-series expansion of 
(Vp" 1 )- That is, the average over $ of d r - is given by 



(d r p- 1 ) = J2[A ± L+B} + 0(p), 



(50) 
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r /M 



6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 



A+ 



-1.964185503296099 
-1.542712134731597 
-1.235264711003273 
-1.008020470281125 
-8.366600265340854 
-7.047957565474786 
-6.014065304058753 
-5.189692421934956 
-4.522475818510165 
-3.975231959999661 
-3.521046167445508 
-3.140087242121197 
-2.817502867825028 
-2.542002591363557 
-2.304886114323209 
-1.500933043143495 
-1.054092553389464 
-7.805574591571754 
-6.011057519272318 
-4.770823620144716 
-3.878143885933070 
-3.214363205318354 
-2.707442873558057 
-2.311598761671936 
-1.996605675599292 
-1.741859372645818 
-1.532923192995987 
-1.359438646187131 



x 10" 2 

x 10~ 2 

x 10" 2 

x 10" 2 

x 10~ 3 

x 10~ 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10" 4 

x 10" 4 

x 10" 4 

x 10" 4 

x 10" 4 

x 10" 4 

x 10" 4 

x 10" 4 

x lO" 4 

x 10" 4 

x 10" 4 



B 



-9.719920769918032 
-7.595781032643107 
-6.072295959309139 
-4.954081856693618 
-4.113353788131433 
-3.467126055149815 
-2.960554843842139 
-2.556541533529994 
-2.229391187912286 
-1.960906506358998 
-1.737934146698723 
-1.550788700414580 
-1.392217662603554 
-1.256707170227638 
-1.140007036978271 
-7.437542878990537 
-5.230319186714355 
-3.876932093890911 
-2.987922074634742 
-2.372891353438443 
-1.929854912525034 
-1.600201924066196 
-1.348310298624353 
-1.151520162987173 
-9.948607781191319 
-8.681200066822863 
-7.641384666121239 
-6.777766234583098 



x 10" 3 

x 10" 3 

x 10" 3 

x 10" 3 

x 10 -3 

x 10 -3 

x 10" 3 

x 10 -3 

x 10 -3 

x 10 -3 

x 10" 3 

x 10 -3 

x 10" 3 

x 10" 3 

x 10 -3 

x 10 -4 

x 10" 4 

x 10" 4 

x 10 -4 

x 10 -4 

x 10" 4 

x 10 -4 

x 10 -4 

x 10" 4 

x 10 -5 

x 10" 5 

x 10 -5 

x 10 -5 



red 



6.8 x 10~ 14 

3.7 x 10" 14 
< 10~ 16 

1.8 x 10~ 14 
1.2 x 10" 14 

8.9 x 10~ 15 
2.5 x 10~ 14 

-2.4 x 10" 14 

8.4 x 10~ 14 
8.9 x 10~ 15 
5.8 x 10" 14 

-4.6 x 10~ 15 

3.5 x 10~ 14 
-3.1 x 10 -15 
-5.7 x 10" 16 
-1.6 x 10" 15 

4.2 x 10" 15 

2.6 x 10~ 14 
3.8 x 10" 15 

9.3 x 10" 15 

1.8 x 10~ 15 

6.9 x 10" 15 

4.2 x 10" 15 
2.0 x 10~ 16 

-2.0 x 10" 16 

1.3 x 10" 15 

2.7 x 10~ 15 
-3.0 x 10~ 16 



-5.0 x 10" 11 
-3.0 x 10" 11 
-5.7 x 10~ 12 
-3.3 x 10" 12 
-7.3 x 10" 12 
-3.7 x 10~ 12 
-1.1 x 10" 11 
6.5 x 10" 12 
-2.9 x 10 -11 
-3.6 x 10~ 12 
-1.7 x 10~ n 
-6.7 x 10~ 15 
-1.1 x 10" 11 
9.4 x 10 -13 
1.8 x 10 -12 



-3 





x 


10" 


14 


-1 


1 


X 


10" 


12 


-6 


5 


X 


10" 


12 


-8 


5 


X 


10" 


13 


-2 


5 


X 


10" 


12 


-2 


3 


X 


io- 


13 


-1 


7 


X 


10" 


12 


-9 


6 


X 


10" 


13 


-8 





X 


10" 


14 


-5 





X 


10" 


14 


-3 


4 


X 


10" 


13 


-5 


6 


X 


10" 


13 



2.9 x 10~ 13 



TABLE III: The first column shows the radius of the orbiting particle in terms of background Schwarzschild coordinate r. 
The second and the third columns show the leading and the sub-leading regularization parameters that we get by a numerical 
matching. The fourth and the fifth columns show the fractional difference between the numerical and the analytic values. 



with and B given by [§] 



.4 



analytic 



B 



(l-SM/ro) 1 / 2 



(1 _ _ 24!) 



analytic 



1/2 



F X 12 

1/2 2 1 



1 1 - 3M/r 
2AI/r 



F 



3/2 



where 



2 r /2 ( M 2 



(51a) 
(51b) 

(52) 



Table III shows the agreement between the regularization coefficients that we get in a radiation gauge and the ones 
in the Lorenz gauge. 



The renormalized self- force is given by f r — ma r , with 

oo 

a ro„r = ^2(a Tetr [£\-A ± L-B) 



(53) 
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ro/M 


TTir ; 2 

F In 


6 


0.031741815 


7 


0.024314669 


8 


0.019541501 


9 


0.016137918 


10 


0.013580536 


11 


0.011595880 


12 


0.010019806 


13 


0.0087455255 


14 


0.0076999148 


15 


0.0068310918 


16 


0.0061012423 


18 


0.0049526422 


20 


0.0040997900 


25 


0.0027292140 


30 


0.0019459393 


35 


0.0014569286 


40 


0.0011313990 


45 


0.00090387369 


50 


0.00073864055 


60 


0.00051979901 


70 


0.00038553381 


80 


0.00029728330 


90 


0.00023619526 



TABLE IV: The radial self-force for a particle in circular orbit around a Schwarzschild black hole. The first column lists the 
Schwarzschild radial coordinate of the orbit; the second gives the radial component of the self-force per unit mass square. 

We can add parts of the perturbation corresponding to a change in mass, and angular momentum in any convenient 
gauge. They, together with an even-parity I = 1 gauge transformation outside tq that accounts for a change in the 
center of mass, have been computed in a Lorenz gauge by Poisson and Detweiler, and we use their results. Note that 
the singular part of the self-force is given by the coefficients in the large- L expansion of the metric and is therefore 
independent of the choice of gauge for I = and I = 1. Our results for the self- force, expressed as the self-acceleration 
a r = fr/xn, are tabulated in Table IIVI 
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IV. GAUGE-INVARIANT QUANTITIES 

The perturbed metric of a particle in circular orbit in a Schwarzschild spacetime is helically symmetric: The helical 
Killing vector k a = t a + il(f> a of the background spacetime Lie derives the particle's trajectory and is a Killing vector 
of the linearized metric /i™J in a Lorenz, Regge- Wheeler, or radiation gauge. As Detweiler noted 0, E3|, the quantity 
H := h a pu a u^ /2 is invariant under gauge transformations generated by a gauge vector £ a that shares this symmetry, 

£ k r = 0. (54) 

A quick covariant derivation follows from the background geodesic equation in the form = and the relation 

u a = u l k a , where u* is the scalar u a V a t. The change in H under a gauge transformation is given by 5^H = u a u"V a £p. 
Then 

uVVo^ = -t»V&,V>Jfe Q = -u^ a Vpu a = 0, (55) 

where Eq. is used in the first equality and the geodesic equation, together with fc a V a (it*) = 0, is used to obtain 
the second. 

In this section we describe the computation of the rcnormalized invariant H and, in effect, compare its values at 
different orbital radii to the values computed in two different gauges by Barack and Sago and by Detweiler. As we 
show in Appendix A of I, the transformation from a Lorenz gauge to the partial radiation gauge we use is generated by 
a helically symmetric gauge vector; the value of H ren in our radiation gauge must therefore coincide up to numerical 
error with its value for a Lorenz gauge. The £ = part of the gauge transformation from a Lorenz gauge to the 
Regge- Wheeler gauge, however, is not helically symmetric, and we must take account of the gauge change in H to 
compare our value to the Regge- Wheeler value. 

As we discuss below, instead of H itself, Sago et al. tabulate a related gauge-invariant quantity that they term 
AU, which is given in terms of H and the background geometry. To facilitate the comparison to their work, we use 
AU as the quantity to tabulate. 

To compute iJ rcn at the position of the particle, we use the harmonic decomposition of the metric to write 

H iet = J2 H imY im (TT/2,0). (56) 

im 

We then match the sequence of values, 

Hl ot = J2 He, n Y em (vr/2, 0) , (57) 

771 

to a power series in £ of the form 

E2 

E ° + (l- 1/2)^ + 3/2) + (I - 3/2)(l ~ l/2)(l + 3/2){£ + 5/2) + " ' ' (58) 

Because the series is obtained from H rct , the renormalization coefficients Ek are again invariant under helically 
symmetric gauge transformations, and we compare in Table [V] the leading term Eq to an analytic form derived by 
Detweiler, 



/ r-3M /ll M \ 

^analytic - ^ r 2 {r _ 2M) ^ g ' ^ 7^2M ) ' ^ 

When H rct is decomposed as a mode sum, the information of the singular part is stored in higher Is. We obtain the 
leading term by matching H\ at to a power series in I from £ ~ 15 to 85. The method of matching is explained in 
Sec. |V] below. 

Table [V] shows the fractional error with which the leading coefficient Eq differs from its analytic form for a set of 
different radii. 



In comparing quantities that we compute in a radiation gauge to those computed by Barack and Sago in a Lorenz 
gauge and by Detweiler in a Regge- Wheeler gauge, we follow the terminology in Sago et al. [if, using the abbreviations 
BS and SD to refer to quantities computed in the two different gauges. We have already mentioned that the comparison 
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ro/M 


\AEo/E \ 


6 


3.8 


X 


io-" 


7 


1.6 


X 


10 -14 


8 


9.0 


X 


10 -16 


9 


9.0 


X 


10 -is 


10 


1.1 


X 


10~ 15 


11 


5.6 


X 


io- 15 


12 


6.0 


X 


io- 16 


13 


2.4 


X 


io- 15 


14 


1.1 


X 


io- 15 


15 


3.4 


X 


10~ 15 


16 


2.7 


X 


io- 15 


18 


4.0 


X 


io- 16 


20 


2.0 


X 


10~ 15 


25 


2.3 


X 


10 -15 


30 


1.4 


X 


IO" 15 


40 


7.0 


X 


10 -16 


50 


1.2 


X 


IO" 15 


60 


1.7 


X 


10~ 15 


70 


< 1(T 16 


80 


5.0 x 10~ 16 



TABLE V: The first column shows the radius of the orbiting particle in terms of the background Schwarzschild coordinate. 
The second column shows the fractional error in Eq, AEq/Eq = \Eo numerical — -Bo analytic] /Eq analytic- 



requires a correction arising from the lack of helical symmetry in the I = gauge vector, but there is an additional 
difference between the BS and SD computations: BS parametrize the perturbed trajectory by proper time r with 
respect to the background metric, while SD uses proper time r with respect to the renormalized perturbed metric. 1 
From the relations 

dr 

u a = u a — , (g aP + h^)u a ul 3 = 1 = g a pu a u^ (60) 
we have, to linear order in the perturbation, 

fir 1 

- = 1 - -h^u a uf. (61) 

Because we, like BS, use proper time with respect to the background metric, we have chosen a Lorenz gauge for 
the 1 = and I = 1 parts of the perturbation so that our u* will coincide with that of BS. In comparing to SD, wc 
then need both corrections - from the I = gauge transformation and from the reparametrization of the trajectory, 
as given in Sago et al. Q. The £ = gauge transformation has gauge vector £ Q = att a , with 

a = m/y/ra(r -3M). (62) 

The relation between u in our modified radiation gauge and in a Lorenz gauge (normalized by proper time with 
respect to the background metric) and u l in a Regge- Wheeler gauge (normalized by proper time with respect to the 
perturbed metric) is then 

u* = u*(l + a - H iea ) + 0(m 2 ), (63) 



1 To maintain a notation consistent with the EMRI literature and with our previous papers, we denote by u a the four- velocity normalized 
with respect to the background metric, with r the corresponding background proper time. Sago et al. use u a and f for these quantities 
and their quantities with no tilde correspond to proper time with respect to the perturbed metric. 
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where 7J ren , given by 



H ien = ^h^u a u p , (64) 



is computed in our modified radiation gauge. Because H lct is invariant under helically symmetric gauge transforma- 
tions, H lcn is similarly invariant. 

Following Sago et al., we write U = u and construct from it another quantity that is invariant under helically 
symmetric gauge transformations by expressing the unperturbed U as a function of SI: With 



«o(fi) == (1 - 3M/R)~ 1/2 , R := ( 



M\ 1/3 

TP) > (65) 



the unperturbed values Uq and Oo satisfy 

U Q =U Q (Sl Q ). (66) 

Then AU is defined as the difference 

AU :=U — Uo(Q.) = U-(1- 3M/R)- 1/2 , (67) 

where the angular velocity that appears in Uq(SI) is the angular velocity of the perturbed trajectory. One can express 
AU in terms of H using the changes SSI and SU in SI and U due to a self-force whose only component is radial, 
f r = ma r , namely 

ssi = si-n a = -ra r ° (1 ~ 2 l^ r °^ a r + o(m 2 ), (68) 

SU = U -U Q = U Q (-H Icn - r -l ar ) + o{m 2 ), (69) 

to obtain 

AU = -(1 - 3M/r )- 1/2 F en . (70) 

Note that AU is gauge- invariant in the standard sense that it has the same value when h a p is replaced by h a p + £^g a /3 
(in this case, when £ Q is helically symmetric). Detweiler and Sago et al. [Hd, [To| also use the term "gauge-invariant" 
to refer to the finite quantities U and f2. Their terminology is motivated by the fact that fl and U are physically 
meaningful: In particular, 51 can be measured by an observer at infinity from the periodicity of received signals sent 
at equally spaced proper times from the orbiting particle. 

One must be careful, however, because SSI and SU are not gauge invariant. With the standard definition of the 
change of a scalar on spacetime under a gauge transformation generated by a gauge vector £™, they change by 

SSI -+ SSI + £jf2, SU -> SU + £ ( U. (71) 

Note, in particular, that Sago et al. use symbols S^Sl and 6%U whose meaning is different from the change in SI and U 
under a gauge transformation, while their symbol S^h a p does have the meaning £^g a /3- We discuss the terminology 
and the reason AU is invariant under helically symmetric gauge transformations in the Appendix. 

Because we have used the Lorenz gauge for the £ = and 1=1 parts of the metric perturbation (derived by 
Detweiler and Poisson) , and the gauge transformation from a Lorenz gauge to the Regge- Wheeler gauge of Detweiler 
is not helically symmetric, we need to make the same gauge adjustment made by BS to our AU to compare its the 
value to that of Detweiler, namely 

AU B = (1 - MI/r )- 1/2 ( r °~ 2M a - H™> KG ) . (72) 

\r -3M J 

In the comparison table, Table IVII below, AU is given in the Regge- Wheeler gauge used by Detweiler. The table 
reports results of the computation performed in the three gauges, radiation, Regge- Wheeler, and Lorenz. 
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r /M 


AU 


A[/(from SD) 


AC/ (from BS) 


6 


-0.29602751 


-0.2960275 


-0.296040244 


7 


-0.22084753 


-0.2208475 


-0.220852781 


8 


-0.17771974 


-0.1777197 


-0.177722443 


9 


-0.14936061 


-0.1493606 


-0.149362192 


10 


-0.12912227 


-0.1291222 


-0.129123253 


11 


-0.11387465 


-0.1138747 


-0.113875315 


12 


-0.10193557 


-0.1019355 


-0.101936046 


13 


-0.092313311 


-0.09231331 


-0.092313661 


14 


-0.084381953 


-0.08438195 


-0.084382221 


15 


-0.077725319 


-0.07772532 


-0.077725527 


16 


-0.072055057 


-0.07205505 


-0.072055223 


18 


-0.062901899 


-0.06290189 


-0.062902026 


20 


-0.055827719 


-0.05582771 


-0.055827795 


25 


-0.043599843 


-0.04359984 


-0.043599881 


30 


-0.035778314 


-0.03577831 


-0.035778334 


40 


-0.026339677 


-0.02633967 


-0.026339690 


50 


-0.020844656 


-0.02084465 


-0.020844661 


60 


-0.017247593 


-0.01724759 


-0.017247596 


70 


-0.014709646 


-0.01470964 


-0.014709648 


80 


-0.012822961 


-0.01282296 


-0.012822962 


90 


-0.011365316 


-0.01136531 


-0.011365317 



TABLE VI: In this table we compare our values of AU with BS and SD. The first column shows the radius of the orbiting 
particle in terms of background Schwarzschild coordinate. The second, third and fourth column shows the values of AU 
computed in a radiation gauge, Regge- Wheeler gauge (SD) and Lorenz gauge (BS). 



V. NUMERICAL MATCHING 



As mentioned in Paper I, we find the singular part of the self-force by matching a power series to its numerically 
computed large- L behavior. Explicitly, we match the sequence of values a rct [£] to successive terms in a series of the 
form 

utlt £w (73) 

where the polynomial, 

k k 

P2k(e) = f[(£ - k - \ + i) f[(£ + k + I - j), (74) 
i=l i=l 

satisfies 



oo 

Y — 



i=0 



= 0. (75) 



When one matches the above scries to the numerical series a Tet [£] over a range £ m i n < I < £ ma x, the accuracy with 
which A and B are are obtained depends on the value of fc max chosen - on the number of parameters used in 
the matching. We choose fc max to minimize the difference (^4(fe max +i) — ^fc ma x)/^femax as a function of fc max - For ^ max 
infinite and no numerical error, ^4(fc max +i) — ^U- max would converge to zero as fc max — > oo. For finite £ max only a finite 
number of parameters can be extracted, and we approximate the value of fc max for which A and B are most accurately 
determined by the value k of fc max that minimizes ^4.(A; max +i) — 4 m „. 

We check this approximation by using our knowledge of the analytic form of the leading renormalization parameter 
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FIG. 1: The solid curve is a plot of |-Eo(fc max +i) — -Eofc max \/Eok max vs fc max . The dashed curve is a plot of \Eoanalytic — 
E f,\/ {Eoanaiytic) as a function of fc max . Both are for orbital radius ro = 10A/ 



Table rvTTI shows for a set of orbital radii the minimum value of |-EWfc max +i) — ^ofc mMC | l^0k max 5 its value at /c max = fc; 
that is compared to the value of \E 0ana iy tic — E k m ^\/ E 0ana i ytic , again at fc max = fc. One can infer from the table that, 
by minimizing ^4(fc max +i) ^ ^4fc max , one obtains values of the renormalization parameters to a fractional precision of 10 -13 
or better. Similar accuracy was reported in Paper I for the leading renormalization parameter in the axisymmetric 
part of ipo- 
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ro/M 


minimum of — £ 




^Ofc max 




alytic 


6 


1 


4 x 10" 14 
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TABLE VII: The first column shows the radius of the orbiting particle in terms of background Schwarzschild coordinate r. The 
second column shows the minimum of -E-o(fc max +i) — £t>fc max |/-Eofc max as a function of fc max . Let that fc max be denoted as k. The 
third column shows the quantity \E 0analytic - E ^\/(E 0ana iy Uc ). £ min = 15. 



The accuracy with which the values of Ek are obtained also depends on the values of £ m in and £ max - For fixed 
f max , there is an optimal value of i m i n that minimizes the error in the regularization parameters. For higher £ m i n , 
contributions of the higher-order E2k are too small to extract; for small l m \ n the i? ret [^] (or o retr [£]) will depart from 
its large- L behavior. Tables IVIIII and IIXI show the behavior of the regularization parameters for different choice of 
train (for a fixed £ m ax = 84) at two different radii. Increasing £ max , of course, increases the accuracy of the computation. 
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TABLE VIII: This table shows the first five fractional differences between successive regularization parameters in the singular 
part of 2H for five different values of £min, with fc max = k. The second, third, fourth, fifth and the sixth columns list the fractional 
differences ( | E n j , x — E n %\/E n ^) for n = 0,-2,-4,-6,-8, respectively. The seventh column gives the fractional difference between 

Eoanaiyticai and BqJ- The last column gives k. All values are for orbital radius ro = 15M. 



The number of renormalization parameters needed to attain machine accuracy is small, because the errors in the 
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TABLE IX: Entries are as in Table [yTTTl with r = 10M. 
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and the |reg parameter| < 1. We work to a fractional error of order 10 13 , terminating the sum over k in Eq. ([75]) at 

Figure [2] shows the result of subtracting successive terms in the numerically-determined large-!/ expansion of the 
self-force, using a retr = / rctr /m, for a particle at tq = 10M. 



ln(residual) 




ln(l) 



FIG. 2: This figure shows a plot of log a r e ct r (and subsequent subtractions of the singular terms) vs log£ The topmost curve is 
of r ~ t. The second curve (from the top) is (a^ ct r - A(£+l/2)) ~ P. The third curve shows (a™* r - A(i+ 1/2) -B) ~ f" 2 . 
The fourth, fifth and sixth curves show the subsequent cumulative subtractions of E2/ -Pa^), Ea/Pa{€) and Eq/P§(1) from the 
third. 



19 



Each successive subtraction should reduce the slope of the \og£ vs log/ retr [£] curve by an integer, and the computed 
graphs shows the expected behavior for the first four renormalization coefficients. 

The accuracy with which the higher-order renormalization coefficients can be computed can be used to estimate 
the accuracy with which iJ ren and / renr are computed. In addition, once the E^k have been found, the terms in the 
mode sum can be grouped as a rapidly convergent numerical sum and an analytically known sum. We have 




B\ 



A(l+ 1/2)- B- 



1/2) — B — 



En 



E% 



E- 



2k u 



Q 



2k u 



E 

1=0 L 



(0 

E 2 

Q2(£) 



E 



where the terms Q2k(£) represent any polynomials in t whose first term is £ 2k . 
allow the second sum on the right side of Eq. (|78j) to be computed analytically, 
used, as in our calculation, the second sum vanishes. 
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The polynomials can be chosen to 
When the polynomials V-zk. can be 



VI. DISCUSSION 



The work here gives a first example of the successful use of a (modified) radiation gauge to compute the self-force on 
a particle in circular orbit. In comparing the renormalized perturbed metric component h a ^u a u^ , we obtain agreement 
to high numerical accuracy with previous calculations in other gauges. The work shows that the singular field can 
be identified with the singular field in a Lorenz gauge. We verify the coincidence numerically to high precision. It 
follows analytically from (1) the expression derived in the companion paper [l[ for the gauge vector £ a that relates a 
Lorenz and radiation gauge and (2) the fact that the gauge transformation of the self-force for a particle in circular 
orbit involves no derivatives of the gauge vector. 

In extending the method to circular orbits in Kerr and then to more general orbits, additional subtleties arise. 
Although spheroidal harmonics decouple in the Tcukolsky equation for ipa and ?/>4, the lack of spherical symmetry 
means that angular harmonics with different values of i no longer decouple in the perturbed metric or in the expression 
for the self-force. This significantly changes the way one computes the contribution to the perturbed metric associated 
with the change in the center of mass - a perturbation that is purely 1=1 for a Schwarzschild background. A second 
complication for non-circular orbits, arises from the existence of a region between periastron and apastron, a region 
where the time and angular harmonics of ipo have a nonzero source. Ways to handle each of these complications will 
be discussed in a subsequent paper. 
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Appendix A: Gauge invariants 

As we mentioned in Sec. IIV1 the symbols 5^fl and d^U used in Refs. [H, 0, [HI in describing U and as gauge- 
invariant are not the standard definitions of the change in a scalar function under a gauge transformation. The 
difference in definitions arises because the four- velocity of a single particle does not conform to the framework used 
to describe perturbations of a set of physical fields on spacetime. The change in the trajectory of a single particle 
involves a four-velocity that is defined only on a single trajectory, and the trajectory is different for the perturbed 
and unperturbed spacctimes. 

Note that Detweiler and Sago et al. are careful to distinguish a lack of gauge invariance of 5fl := fl — £Iq and 
5U := U — Uq from the invariance of 51 and U . The number CI can be measured, for example, by an observer at 
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infinity from waves received from the orbiting particle: It is in this sense gauge-invariant. However, when f2 is regarded 
as a function of trajectories or as a function on a Schwarzschild spacetime that assigns to each point of spacetime the 
frequency of a circular orbit through that point, it becomes gauge-dependent. In particular, in the active description of 
a gauge transformation, an infinitesimal diffeo changes the metric at a point of spacetime by £^g a /j, and the frequency 
of a geodesic through that point then changes by £^f2. (In the passive description of a gauge transformation, with 
f2 regarded as a function of the coordinates, an infinitesimal coordinate transformation correspondingly changes the 
functional dependence of Q on the coordinates.) 

We give here a brief review of general relativistic perturbation theory and gauge transformations. We describe 
linearized perturbation in terms of a family of fields, with gauge transformations described in the active sense of 
diffcos (diffcomorphisms) . We begin with a family of metrics <? Q ^(A), defining 5g a p as the change in the metric to first 
order in A as the quantity 



5 9af3 = jTdaflW 



(Al) 



A=0 



When no other fields are present, the metric g a p is physically equivalent to the diffeomorphically related metric i>* g a p, 
where ip is a diffeo (a diffeomorphism) and ip* is the pullback, whose components are given by 



dib a dib T 



(A2) 



A gauge transformation is defined by considering a smooth family of diffcos with ipo the identity. We denote 
by the vector field tangent at each point P to the orbit A n- ip\{P) of P. That is, to linear order in A, a point 
with coordinates x^ is mapped to a point with coordinates ip^(P) = x^ + A£ M (P). The family of metrics g a p{X) is 
then physically equivalent to the family ?/>^g Q ^(A), and the metric perturbation 8g a p is physically equivalent to the 
gauge-related metric perturbation 



a=o 



A=0 



Jaf3 



(A3) 



A=0 



When other physical fields [Ti , . . . , Tjv] are present, the families of fields \g a p(X), Ti(\), ■ . . , T/v(A)] and 
bP\9a/3, V'a-^Ij ■ • ■ ' V'a-^v] describe the same physical system. The corresponding physical equivalence of the perturbed 
fields, linearized about A = 0, is the equivalence of ST and ST + £^T, for each physical field T, where 



ST=-T(X) 



A=0 



= ST + S^T, and <5 € T = £ € T. 



A=0 



In particular, if T is a scalar, ip^T = T o implying 

T(P) = rxT[^\P)]. 
A scalar / constructed locally from a set of tensor fields Ti, . . . Tk, with 



(A4) 



(A5) 



f = F(T 1 ,...T k ), with /(P) = FpUP), . . .T fc (P)], 



(A6) 



satisfies 



4>f = FWr 1 ,...fpT k ), 



(A7) 



implying, for any vector field £ Q , 



A£^Ti, . . . , Tfc + A£^Tfe) 



A=0 



(A8) 



One also uses a set of reference fields - basis vectors or coordinates, for example - that are independent of A and 
are not changed by the diffeo that maps the set of physical fields to the physically equivalent set. For example, with 
t a the vector 9 t , the gauge transformation S^gtt of a component g a pt a t^ is 



(A9) 



We now turn to the problem at hand, the behavior of U and Vt under gauge transformations. These are defined 
in terms of the four- velocity u a . To define them and u a as fields, to write a change Su a as a vector at a point using 
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the definition (|A4I) . and to define the Lie derivative £^u a , one must introduce a set of nearby orbits. In the problem 
considered here, one can define an unperturbed vector field u a by taking u a (P) to be the four- velocity of the circular 
orbit through P. The perturbed vector field at P is then the four- velocity of the perturbed circular orbit through P. 
The scalar U of Sec. IIVI is constructed locally from g a/ 3,H,,t a , and 4> a : 

u = u[ 9aP (t a + n<j> a )(te + n^)]- 1 ' 2 =■. u^si, t a , <p a ). (aio) 

Here t a and (j) a are reference fields, defined as the vectors dt and d<p tangent to the coordinate lines, so that they have 
the same values for the perturbed and unperturbed metric. (Note that they are not Killing fields of the perturbed 
metric, and cannot be defined in this way.) If is helically symmetric with respect to k a = t a + VL(P)(f> a , for a given 
fixed P, then 



hU{P) = ^U{g al} + \£ i9afi ^ + \£^,t a ,4> a ) 



= £ € U{P). (All) 

A=0 



Similarly, 

6^n(P) = £ € 0(P). (A12) 

Finally, defining Uq as in Eq. (|66[) . we have 

S € [U-U (n)]= [£ € U-£ € U o ]\ x=o = £ € U o -£ € U o = 0. (A13) 

That, is, U — Uq(Q) is gauge invariant in the usual sense, for helically symmetric gauge vectors. 

If we denote by 5^U and the changes in U and introduced by Sago et al., Eqs. (|All[) and (|A12[) are equivalent 
to writing 

S^U = 0, 5 € n = 0. (A14) 

Although any scalar constructed locally from physical fields alone has this behavior, Eqs. (|A11[) and (|A12[) (or, 
equivalently, (|A14[0 are not trivial relations, because they involve reference fields that do not change under gauge 
transformations. 

The difference between 5^ and 5^ is the change associated with a second kind of gauge freedom, mentioned in 
Appendix A. 2 of Paper I, namely an infinitesimal change (by a displacement £ Q ) in the background geodesic to which 
one compares a geodesic in the perturbed spacetime. 
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